Noise-seeded spatio-temporal modulation instability in normal dispersion 
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In optical second harmonic generation with normal dispersion, the virtually infinite bandwidth 
of the unbounded, hyperbolic, modulational instability leads to quenching of spatial multi-soliton 
formation and to the occurrence of a catastrophic spatio-temporal break-up when an extended beam 
is let to interact with an extremely weak external noise with coherence time much shorter than that 
of the pump. 

PACS numbers: 



The noise-seeded instability of extended wave pack- 
ets (WP) in conservative evolutional non-linear systems 
is a general and relevant phenomenon in wave physics, 
whose main manifestations are the appearance of regu- 
lar modulations, wave break up and, eventually, local- 
ized or soliton-like substructures. The theory used for 
the description and interpretation of the resulting rich 
phenomenology is that of the modulational instability 
(MI) of plane and monochromatic waves, extensively in- 
troduced in the context of gravity waves in deep waters 
and applied for several different systems including 
plasmas |2j, electric circuits 1311 , Bose-Einstein conden- 
sate y and, of course, optics |5|]. The usually addressed 
MI signature is the preferential noise amplification at 
a given, intensity dependent, (spatial or temporal) fre- 
quency, which causes regular modulation in the direct 
space and side-bands in the spectral domain. This fea- 
ture suitably describes the MI of both mono-dimensional 
(ID) systems and of multi-dimensional "elliptical" ones, 
i.e. those supporting equi-sign linear phase modulation 
in all the available dimensions. However, it is gener- 
ally not adequate for "hyperbolic" systems, where op- 
posite signs occur for different dimensions. The ellip- 
tical is the most frequently encountered regime in case 
of matter waves in isotropic media. The hyperbolic, in 
contrast, is the typical case of optical WPs in normally 
dispersive bulk media, diffraction and chromatic disper- 
sion leading in this case to linear phase modulations with 
opposite signs. Recently, dispersion-management tech- 
niques based on use of periodic potentials |6j have made 
the hyperbolic regime of great interest also for the Bose- 
Einstein Condensate (BEC) waves. 

As clearly pointed out in the first theoretical analysis 
of hyperbolic MI performed by Luther et Al. for Kerr 
non-linearity in optics, and also evident from the anal- 
ysis of the X( 2 )-driven MI ||, the the key feature that 
distinguishes the hyperbolic (normal) from the elliptic 
(anomalous) instability regime is that the MI gain profile 
in the k-to space is unbounded in the first (and only in the 



first) case. Indeed, in the frame of the usually adopted 
parabolic approximation for the material dispersion, we 
should say that any fluctuation with arbitrarily large spa- 
tial and temporal frequency shift respect to the carrier 
mode has to be amplified, provided that both shifts lye on 
the suitable hyperbolic surface in the k-u: domain. 

The unbounded feature of the hyperbolic MI rises two 
relevant questions: (i) the first one concerns the interpre- 
tation of the number of studies reguarding spatial MI in 
multi-dimensional systems. Indeed, following the first ex- 
perimental MI demonstration in (anomalous dispersive) 
ID-temporal optical fibers a number of experiments 
have been performed in multi-dimensional (i.e. in planar 
wave guides or in bulk samples) normal-dispersion ma- 
terials, addressing the ID spatial break up of extended 
beams driven by quadratic {lflL Ull ] as well as cubic 
ultrafast non-linear response. Surprisingly enough, the 
results were successfully interpreted in terms of the di- 
rect spatial analogous of the temporal MI of above, as 
if the temporal degree of freedom and so the hyperbolic 
nature of the instability had not taken any part in the 
process, (ii) The second question concerns the scenario 
that one should expect when the system interacts with 
a very broad band (i.e. virtually (5-correlated) noise. In 
this case, in fact, the unbounded feature of the MI gain 
should lead one to forecast a catastrophic break in the 
ST domain (e.g. down to the numerical grid in calcu- 
lations, in the quoted approximation), no matter how 
weak the input noise is. We note that the possible oc- 
currence of such catastrophic dynamics has never been 
considered in the literature. In fact, hyperbolic MI has 
been studied only for the case of bell-shaped, noise-free, 
input wave packets, e.g. for investigating the impact of 
ST self-phase-modulation (SPM) on filament formation, 
pulse splitting and related phenomena 0. To the best 
of our knowledge, the role of the noise has been con- 
sidered only in the context of ID models. The aim of 
this work is that of providing experimental evidence of 
the genuine hyperbolic feature of the noise-seeded MI in 




FIG. 1: Calculated space-time profiles (bottom figure) and corresponding integrated fluence profiles (top figure) for the fun- 
damental harmonic propagated in a 50mm lithium triborate (LBO) crystal in regime of Ak — 2k(uo) — k(2uio) = 5cm -1 , in 
absence of noise (a and d), and with two different kinds of noise (b and e): 6nm and 25mrad of bandwidth, (c and f): 20nm 
and 50mrad. The intensity level of noise is 1%. The input is a gaussian FH WP with lps FWHM duration, 600^im FWHM 
beam width, and 40GW/cm 2 peak intensity. Calculations are performed in the frame of ID (spatial) + ID (temporal) + ID 
(propagation) model. Asymmetry in the temporal coordinate is due to group-velocity mismatch. 



bulk, normally dispersive, optical systems. To this end 
we performed experiments and calculations in which, for 
the first time to our knowledge, a controlled, broad-band 
noise is injected together with the strong, (quasi) plane 
and monochromatic pump into the system. 

The particular system that we have chosen to investi- 
gate is the same as in Ref. |l3|. i.e. that of an optical 
WP shaped as a large (with respect to diffraction) , elon- 
gated beam and long (with respect to dispersion) pulse 
that propagates in a non-linear crystal tuned for 

second harmonic (SH) generation close to phase match- 
ing. For the chosen crystal (lithium triborate, LBO) and 
wavelength (first harmonic, FH, 1055nm) the chromatic 
dispersion (see caption to Fig. 1 for details) is such that 
one should expect MI to take place in the hyperbolic 
regime Q. We performed numerical calculations in the 
frame of 2D+1 model integrating (via FFT, split step 
and Runge Kutta algorithms, with up to 15fs 6/im grid 
and 40/xm step) the x* 2 ' coupled-wave equations for FH 
and SH envelopes Ej(z, x, t) 

L(uj ) E 1 + X E 2 El expHAfcz) = 0, . . 

L(2u) Q ) E 2 + i5Vd t E 2 + xE\ exp(iAfcz) = 0, 1 ' 

where L(uj) = id z + (2k(uo)y 1 d 2 xx ~ (k"(u)/2)d? t , k" is 
the group-velocity dispersion (GVD), SV = fc'(2o) ) — 



k'(u>o) weighs the group-velocity mismatch (GVM), and 
Ak = 2fc(cj ) — k(2uj )- Figure 1 gives the calculated 
fluence (e.g. energy density) profile of the FH at the 
output of a 50 mm crystal (top) , together with the corre- 
sponding spatio-temporal (ST) intensity maps (bottom). 
Figs. la,d, which refer to noise-free input, show the oc- 
currence of a regular, highly contrasted, spatial break up 
of the beam into a spatial-soliton array, which appears as 
the consequence of MI seedin g by the deterministic wave- 
envelope modulation (WEM 13]). When the ST noise is 
injected, WEM and noise-seeded MI compete and the re- 
sults are those shown in Fig. lb,c, and Fig.le,f, where 
the input-noise bandwidth (BW) is increased (from left 
to right) while keeping fixed the noise intensity at 1% 
of the level of the pump. See how the impact of the 
noise dramatically increases on enlarging its BW, owing 
to the unbounded instability. Note how the noise, instead 
of deepening the spatial modulation (as it occurs in the 
frame of the 1D+1 models ^(J), quenches it almost com- 
pletely. The reason is the appearance of a "chaotic gas" 
of localized ST structures (Fig. If), which gets washed 
out by the temporal integration. 

In order to verify if the outlined catastrophic behavior 
is a genuine physical effect or an artefact of the approx- 
imations adopted we performed a SH-generation labo- 
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FIG. 2: Measured fluence profiles (top figure) of FH WP recorded by a CCD camera at the output of the 50mm-LBO crystal, 
and (bottom figure) corresponding beam profiles (note that the horizontal scales of images and profiles are not the same) in 
conditions: (a) and (d), with no noise; (b) and (e) with noise intensity 0.01% of that of the pump; (c) and (f), with 0.1% of 
noise. The noise spatial and temporal BWs were respectively of 60mrad and lOOnm. The input pump intensity was 20GW/cm 2 . 



ratory experiment in similar conditions to those which 
Fig. 1 refers to. To this end we used a strongly elon- 
gated (1000x70 fJ,m), long (lps) pump WP, as clean as 
possible from any spatial or temporal substructure, pro- 
vided by a CPA Nd:Glass laser (TWINKLE, Light Con- 
version). Then we superimposed to the pump a weak, 
broad-bandwidth, ST noise, of controllable intensity, gen- 
erated on a separate channel by a broad-band quantum- 
noise parametric amplifier. For the noise generation we 
used a 15mm LBO crystal pumped by the SH of (a por- 
tion of) our pump pulse. Both pump and noise WPs 
were launched synchronously into a 50 mm LBO crystal, 
tuned for phase matched SH generation. The spatial and 
temporal BWs of the noise field at the input of the SH 
generator were 100 nm and 60 rad, respectively. Fig. 2 
(top) shows the fluence distributions of the FH beam at 
the crystal output facet as recorded by a CCD camera 
and suitable imaging optics. The corresponding profiles 
(along the long axis of the beam) are given in Fig. 2 
(bottom) for a more quantitative description. The re- 
sults in the left, center and right part of the figure refer 
to average noise fluence 0, 0.01% and 0.1% of that of the 
pump, respectively, for a fixed noise ST BW. The result- 
ing scenario fully confirms the model prediction. Indeed 
the noise-induced quenching of the WEM-seeded spatial 
MI takes place in the experiment for a lower noise level 
than in calculations, which indicates that the accessible 
BW is even larger than the computational BW used for 
obtaining the Fig. lc,f results. 

The described, near-field, measurements have the obvi- 
ous limitation of confirming the model prediction only on 
the basis of a time-integrated effect, the underlying ST 
structure not being detectable by any technique. More- 
over, both the numerical and the experimental results 



that we have presented do not produce a direct evidence 
of the hyperbolic nature of instability. In what follows, 
in order to overcome these limitations, we illustrate the 
complementary, far- field analysis, concerning the charac- 
terization of the angular spectra (AS) of the generated 
field. Figure 3a contains the calculated AS (i.e. the 
square modulus of the field Fourier transform) of a FH 
profile analogous to that in Fig. If (see the caption for 
details). The bright central spot corresponds to the spec- 
trum of the input pulse while the surrounding structure 
describes the amplified fluctuations. Note the evident hy- 
perbolic shape of the instability region, which coincides 
with the region where the calculated Mi-gain profile is 
the largest |g. We verified that, no matter the size of the 
Fourier space, the instability always reaches the border of 
the spectral box (in case of large enough noise BW), thus 
confirming the unbounded nature of the MI process in 
the frame of the adopted model. For the far-field exper- 
iment, we measured the FH AS by placing the entrance 
slit of a large-numerical-aperture imaging spectrometer 
at the focal plane of a positive lens. The detected AS in 
the short-wavelength branch of the spectrum (the long- 
wavelength one is not accessible due to sensitivity cut off 
of the silicon CCD detector) is reported in Fig. 3c. The 
detected portion of the AS (the radiation at large angles 
was clipped by the very narrow aperture of our 3x3x50 
mm 3 LBO crystal) exhibits a good qualitative agreement 
with calculations (see the zoom in fig. 3b) and confirms 
the genuine hyperbolic feature of the instability process. 

In conclusions, in regime of second- harmonic gener- 
ation with normal dispersion, by superimposing on a 
intense, clean, pump wave packet and a weak spatio- 
temporal noise (with coherence time much shorter than 
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FIG. 3: a) Calculated angular spectra (AS) of the FH field 
for operating conditions as those in Fig. If, but for larger 
input-noise BW (namely: 40nm and 50mrad); (b) Zoom of 
(a) corresponding to the region of detection; (c) measured AS 
of the radiation exiting the crystal; (d) measured AS of the 
input noise. Grey colors refer to logarithmic scale. 



that of the pump) we have shown that the noise-seeded 
MI develops in the spatio-temporal domain. Due to the 
unbounded feature of the instability, and so the large re- 
sponse of the system to white noise, a noise as weak as 
0.01% of the pump leads to quenching of the spatial de- 
terministic beam break up and spatial-soliton formation 
caused by wave envelope modulation (WEM). The analy- 
sis of the results in the spectral domain outlined the gen- 
uine hyperbolic feature of the instability, which couples 
different frequencies to different angles. Because of this 
coupling, the quenching of detectable spatial effects can- 
not be simply interpreted as the averaging of several in- 
dependent spatial structures, occurring for different time 
slices of the wave packet. As pointed out in ^jji the typ- 
ical (amplitude and phase) fluctuations that were trig- 
gering the instability in previous experiments (without 
external noise injection) were probably caused by laser- 
beam, or optical-component or non-linear sample imper- 
fections, thus leading to a "frozen-noise" with the same 
coherence time of the pump (that virtually coincides with 
the pulse duration, all lasers operating close to the trans- 
form limit). This might explain why the resulting in- 
stability was so successfully described in the frame of 
monochromatic models. Finally, we expect that the de- 
scribed, hyperbolic instability should play a dramatic role 
when the field is strong enough to probe vacuum state 



fluctuations, which indeed provide the source of a virtu- 
ally (^-correlated noise. We expect that quantum-noise 
seeded hyperbolic instability should dominate not only 
the X^ 2 ) parametric-amplification regime (as evident, for 
example, in ref. |l4| and also in the more recent rcf. 
(lif) but also the classical, unseeded, second-harmonic 
generation process and, possibly, the Kerr regime too. 
This might explain the " spontaneous" quenching of spa- 
tial break-up seen at high pumping in 0, E5 • Owing 
to the unbounded nature of the instability, the robust- 
ness of non-linear dynamics of normally dispersive media 
with respect to the interaction with the quantum noise 
represents therefore a crucial issue with deserves further 
investigation. 
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